Bilayer graphene Origami: curvature-induced p-n junctions 
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A massive quantum particle on a two-dimensional curved surface experiences a surface-geometry 
induced attractive potential that is characterized by the radii of curvature at a given point. With 
bilayer graphene sheets and carbon nano-ribbons in mind, we obtain the geometric potential Vg 
for several surface shapes. Under appropriate conditions that we discuss in detail, this potential 
suppresses the local Fermi energy. Therefore, we argue that in zero band-gap materials with a 
quadratic band structure, it will create p- and n-type regions. We discuss the consequences of this 
result, and suggest that surface curvature can provide a novel avenue to create p-n junctions and, in 
general, to control local electronic properties in carbon nano-ribbons and bilayer graphene sheets. 

PACS numbers: 72.80.Rj,02.40.-k,71.10.Pm 
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Introduction: The classical dynamics of a particle sub- 
ject to (holonomic or non-holonomic) constraints is well- 
understood. In Lagrangian formulation, these con- 
straints can be represented by an effective potential^ In 
quantum theory, the potential representing constraints 
plays an instrumental role in determining the quantized 
spectrum of the system* 2 " Indeed, low-dimensional sys- 
tems are realized by constraints that classically restrict 
the motion of the particle in one or more directions. For 
example, in two-dimensional electron (or hole) systems 
(2DES), typically realized in ffl-V semiconductor het- 
erojunctions^ the barrier potential at the hetero junction 
interface, created by differing electron affinities, leads to 
a clear separation of energy scales for the motion in the 
2D plane and motion along the direction normal to the 
plane. It allows one to separate the particle wavefunction 
into normal and surface components and, at low energies, 
focus solely on the surface component A 

The quantum properties of a non-relativistic particle 
constrained 5 on an arbitrary orientable surface were first 
investigated by da Costal He showed that apart from 
the necessity of expressing the Laplace-Beltrami opera- 
tor V 2 in terms of the curvilinear surface co-ordinates, 
the geometry of the surface also leads to an attractive 
geometric potential Vg(qi, Q2) — —{h 2 /8m)(Ki — k 2 ) 2 
where m is the mass of the particle, (qi, q%) denote surface 
co-ordinates, and Ki,kq are the two position-dependent 
principal curvatures of the surface. Note that this poten- 
tial is purely a result of particle confinement, and there- 
fore is independent of the electric charge of the particle; 
it is therefore the same for electrons and holes. This re- 
sult is applicable in the limit Wn — > where W is the 
thickness of the "two-dimensional" surface and k is the 
surface curvature. W will correspond to the width of the 
quantum well in 2DES. The geometric potential is small 
near a surface maximum or minimum, and is large near 
saddle regions where K1K2 < 0. In particular, we note 
that for surfaces with a constant curvature (spherical sec- 
tions or a plane) the geometric potential is equal to zero, 



Vg = 0. The bound states that appear due to the geo- 
metric potential® and the resultant resonance microwave 
absorption 9 have been investigated in twisted^ or bent 
(quasi one-dimensional electron) waveguides^ The effect 
of curvature on spin-orbit coupling has been studied in 
the case of nanospheres and nanotubesii 2 - More recently, 
the geometric-potential induced charge separation in he- 
licoidal ribbons has been analytically explored. 13 How- 
ever, absent a truly two dimensional system that can be 
(easily) bent, the effect of geometric potential on the elec- 
tronic band structure has been justifiably ignored. 

In this paper, we point out that the geometric po- 
tential provides a hitherto unexplored avenue to locally 
change the carrier polarity in zero-gap materials. We 
predict that it will be possible to create p-n or n-p-n 
junctions purely from a suitable surface geometry. Bi- 
layer graphene provides an ideal realization of a zero- 
gap material with chiral electrons and holes that have 
identical mass. Graphene nano-ribbons (GNR), recently 
created by unzipping single and multi-walled carbon nan- 
otubesi^ may provide another promising candidate that 
can be used to explore the curvature effects in a gapless 
semiconductor. Therefore, in this paper we primarily fo- 
cus on bilayer graphene. In the next section, we obtain 
the geometric potential for two generic cases: first a sur- 
face z = f(x) that is isometric to a plane and second 
a catenoid that is not isometric to a plane. In the fol- 
lowing section, we show that at present carrier densities 
n2D, carriers in bilayer graphene, with its Bernal-stacked 
hexagonal lattice structure, can be treated as massive 
particles on a curved surface. We conclude the paper 
with a qualitative discussion about experimental conse- 
quences and speculation regarding the geometric poten- 
tial for monolayer graphene with massless chiral carriers. 

Monolayer and bilayer graphene have been studied in 
great detail in recent years. P-N junctions in monolayer 
graphene are fabricated by local top-gating^ and with 
an air-bridge top gate^ 6 ' Bilayer graphene has been ex- 
plored for its electric field tunability i 17 i 18 The effect of 
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random surface curvature (ripples) on electronic prop- 
erties in monolayer graphene is primarily described via 
gauge potentials^; however, the effect of geometric po- 
tential in monolayer and bilayer graphene has not been 
investigated. In monolayer graphene near the neutral- 
ity point, disorder induced electron-hole puddles and the 
resulting array of p-n junctions is observed.— However, 
in this case, the origin of the disorder and the role, if 
any, played by monolayer surface curvature is unclear^ 
In particular, the topic of p-n junctions made with bi- 
layer graphene (that we propose in this paper) has been 
scarcely investigated. 

In this paper, we argue that due to its essentially 
two-dimensional nature (Wk — ► 0), gapless quadratic 
band structure, and exceptional material strength, bi- 
layer graphene presents an excellent candidate for explo- 
ration of two-dimensional systems with curvature. 
Continuum Model: Let us consider a non-relativistic par- 
ticle of (effective) mass m constrained to a surface given 
by z = f(x) (See Fig. [I}. We characterize a point 
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FIG. 1: (Color online) Schematic shapes of various curved 
surfaces considered in this paper. The top two, a sheet with 
a cylindrical bump with radius a (a) and a sheet suspended 
between two mesas (b), have only one nonzero principal cur- 
vature k{x) and are isometric to a plane. The surface at the 
bottom, a catenoid (c), has two nonzero principal curvatures 
at each point Ki(x,y)K2(x,y) < 0. 

on the two-dimensional surface embedded in the three- 
dimensional space by r = (x,y, f(xj). The geometric 
potential is obtained from the two fundamental forms 
associated with the surfaced The first fundamental form 
or, equivalently, the (diagonal) metric tensor— for the 
surface is given by = d^r ■ djV — dia(l + f' 2 (x),l) 
where i,j — x,y and /' = df/dx. The second funda- 
mental form hij = —diT ■ djii characterizes the change 
in the surface-normal n along the surfaced and has only 
one nonzero element, h xx {x) = —f"(x)/[l + f'(x) 2 } 3 ^ 2 
where /" = d 2 f/dx 2 . Principal curvatures of the surface 
are eigenvalues of the matrix^ a = —hg^ 1 and, in this 
case, they are given by K\(x) = f"(x)/y/i + f'(x) 2 and 
k-2 = 0. The zero eigenvalue reflects the translational 
invariance of the system in the y-direction, and the zero 
Gaussian curvature K — n\K2 = everywhere on the 



surface implies that it is isometric to a plane. 21 Thus, 
for such a surface, the geometric potential is given by 
V G {x) = -{H 2 /8m)[f"(x)/y/l + f'(x) 2 } 2 . _ 

We remind the reader that the potential is purely quan- 
tum, can be different for isometric surfaces (for example, 
a plane and a cylinder), and is always attractive. The 
Hamiltonian for the constrained particle is H = Hq + Vg 
and the kinetic contribution is given by 
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where g = det <jty, and s(x) — J* du^Jl + f'(u) 2 is the 
arc-length along the surface. In terms of the surface co- 
ordinates (s,y), the geometric potential term becomes 



Vc(s(x),y) 



8m 



s"(x) 



s'{x) 2 ^/s'(x) 2 - 1 



(3) 



Note that for a cylindrical bump of radius a (Fig.QJi), we 
recover— a constant geometric potential Vg = —h 2 /8ma 2 . 
This attractive geometric potential is independent of 
the sign of the curvature: a bump, a trough, or com- 
binations thereof lead to the same constant attractive 
potential. For a suspended sheet (Fig. QJ>)j 22 f(x) — 
a [1 — cosh(x / a)] where a sets the length-scale for the 
catenary, the geometric potential is given by 



V G (s(x),y) 



8ma 2 



(4) 



Lastly, we consider the particle confined to a catenoid 
(Fig. [TJ;) parameterized by two dimcnsionlcss surface 
co-ordinates {u, v) where x = acosh(u) cos(u), y — 
acosh(u) sin(w) and z = au. Note that perfect (bilayer) 
graphene sheets are planar and, as such, they cannot be 
mapped onto a surface with nonzero Gaussian curvature. 
However, fullerenes or dislocations in the honeycomb lat- 
tice can lead to surfaces with nonzero Gaussian curva- 
ture. A straightforward calculation gives the geometric 
potential 



V G (u,v) 



1 



2ma 2 cosh 4 (w) 



(5) 



We note that Vg(u,v) is solely a function of u, just as 
V G (s,y) in Eq.Q is solely a function of s. However, 
unlike in the former case, the kinetic energy term on the 
catenoid H = — (h 2 m(u) /2ma)(d 2 / du 2 + d 2 / dv 2 ) is not 
separable and reflects the nonzero Gaussian curvature of 
the surface. 

Applicability to Bilayer Graphene: To map a curved bi- 
layer graphene sheet onto a quantum particle on a curved 
surface, two criteria must be satisfied: kpaQ <C 1 and 
naG 1 where a G — 1.4 A is the carbon-carbon dis- 
tance in graphene, k,F is the Fermi wavevector, and k 
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is the surface curvature. The first criterion ensures that 
bilayer graphene can be treated in the continuum limit. 
In this limit, for carrier densities ri2D ^ 3.5 x 10 13 cm" 2 , 
only the gapless valence and conduction bands are occu- 
pie d 23 ' 24 and the carriers can be approximated as parti- 
cles with an effective mass m* ~ 0.03m e . The second 
criterion ensures the detailed lattice structure can be 
ignored on the curvature length-scale, and because the 
"thickness" of bilayer graphene sheets is W ~ 3a G , it 
also ensures that Wk <C 1, or that the bilayer graphene 
sheet with nonzero thickness can be considered a sur- 
face. The first criterion is easily satisfied at present car- 
rier densities in bilayer graphene with n2D ~ 10 12 cm~ 2 
and corresponding kp ~ 10~ 2 A _1 . The second crite- 
rion, kclg <C 1 is satisfied by carbon nanotubes with typ- 
ical radii 1/k ~ 10-30 A and by typical surface ripples 
on monolayer graphene.— When k ~ kp the geometric 
potential Vg becomes comparable to the Fermi energy 
Ep = h 2 k F /2m* , whereas in "planar" bilayer graphene 
k w < k F < l/a G . 

Within the continuum approximation, we can consider 
two cases. When l/oG>K>/ijr the problem is reduced 
to that of a single-particle in an attractive external po- 
tential. 13 This case is not very realistic since at (vanish- 
ingly) low carrier densities, disorder effects dominate the 
physics . 20 ' 26 On the other hand, when k ~ kp <C 1/ac, 
the geometric potential will reduce the local Fermi en- 
ergy. In particular, if we consider graphene with n-type 
carriers in the flat region (Epo > 0), the geometric poten- 
tial will lower the Fermi energy from the conduction band 
to the valence band, Ep(qi,q2) = Ep — \ Vg(qi, Q2)\ < 0. 
Note that as long as the geometric potential is constant 
or slowly varying (as would be the case for a cylindri- 
cal bump), this result is robust and is independent of 
any other quantum numbers, such as chirality or elec- 
tric charge, that the carriers may have. Therefore, we 
predict that in such a configuration, the curved regions 
will have p-type carriers with a natural p-n junction cre- 
ated between the curved and flat regions. For the surface 
2 = f(x), the geometric potential will create p- and n- 
type strips where the location, size, and polarity of each 
strip is determined by the local Fermi energy Ep(x). 

Bilayer graphene is distinguished from its semiconduc- 
tor counterparts by the chiral nature of its carriers that 
leads to zero transmission probability for normal inci- 
dence at a potential barrier in spite the zero gap be- 
tween the electron and hole bands^i Thus, the chiral- 
ity of the carriers in bilayer graphene will only affect 
the transparency (i.e. the transmission amplitude as a 
function of incidence angle) of the p-n junctions created 
by the geometric potential* 2 ^ Such p-n junctions can be 
created by cither depositing bilayer graphene on a pat- 
terned substrate with requisite bumps and troughs with 
radii R ~ 100 A or by suspending it across a narrow 
channel. A local probed may also offer an alternative 



way to induce controlled curvature in a suspended bi- 
layer graphene sheet. We emphasize that our prediction 
is valid for any two-dimensional zero-gap semiconductor 
with appropriate curvature, although bilayer graphene 
and GNRs are the most promising candidates. 
Discussion: In this paper we have proposed that bilayer 
graphene on a patterned substrate provides a novel av- 
enue to create a p-n junction. Our prediction provides 
a direct way to experimentally probe the effect of ge- 
ometric potential that has been, hitherto, neglected in 
two-dimensional systems. 

The origin of the geometric potential Vg is well- 
understood for a non-relativistic massive particle. 6 In 
case of a planar 2DES, a lattice model naturally leads 
to carriers with an effective mass. However, a lattice- 
model generalization of the geometric potential is, to our 
knowledge, an open question. The primary difficulty in 
addressing such a question is that typical lattice models 
start with zero thickness, W — 0, whereas the geometric 
potential arises from taking the limit W — •> 0. 

We conclude the paper with speculation about curved 
monolayer graphene. Since the geometric potential Vg 
is obtained from non-relativistic Schrodinger equation, 
prima facie is not applicable to monolayer graphene 
where the carrier dynamics is well-described by the two- 
dimensional Dirac Hamiltonian. Monolayer graphene 
with random curvature (ripples) is modeled using ran- 
dom gauge fields^ 9 - and mean-curvature-dependent po- 
tential.^ - We are unaware of a rigorous derivation of ge- 
ometric potential for massless particles or for a lattice- 
model that leads to linearly dispersing bands. Dimen- 
sional analysis implies that for a surface z — f(x), the 
geometric potential, if nonzero, must scale as V G oc 
7m;g|k| where v G is the velocity of massless carriers in 
graphene.— It follows that under the conditions discussed 
in the last section, this geometric potential will locally 
change the Fermi energy Ep(x) and will lead to the for- 
mation of a p-n junction. We remind the reader that the 
transparency of such a junction in monolayer graphene 
is drastically different from that in bilayer graphene al- 
though both are zero-gap materials^ 

The test of our prediction via experiments and the 
generalization of the geometric potential to a lattice 
model will provide insights into the electronic structure 
of curved two-dimensional materials with zero band-gap. 
Curvature induced p-n and p-n-p junctions in bilayer 
graphene and the GNR may open up new vistas for nano- 
electronics. 
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